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ABSTRACT 

In this paper we take a deeper look at the technically ele- 
mentary but physically robust viewpoint in which the Casimir 
energy in dielectric media is interpreted as the change in 
the total zero point energy of the electromagnetic vacuum 
summed ov er all states. Extending results pr esented in previ- 
ous papers [ ^iep-th/9609T9^ ; ^iep-th/97020& 7f we approximate 
the sum over states by an integral over the density of states 
including fimte volume corrections. For an arbitrarily-shaped 
finite dielectric, the first finite-volume correction to the den- 
sity of states is shown to be proportional to the surface area of 
the dielectric interface and is explicitly evaluated as a function 
of the permeability and permitivity. Since these calculations 
are founded in an elementary and straightforward way on the 
underlying physics of the Casimir efi'ect they serve as an im- 
portant consistency check on field-theoretic calculations. As 
a concrete example we discuss Schwinger's suggestion that 
the Casimir effect might be the underlying physical basis be- 
hind sonoluminescence. The recent controversy concerning 
the relative importance of volume and surface contributions 
is discussed. For sufficiently large bubbles the volume effect 
is always dominant. Furthermore we can explicitly calculate 
the surface area contribution as a function of refractive index. 

PACS: 12.20.Ds; 77.22. Ch; 78.60. Mq 



I. INTRODUCTION 

The Casimir effect in dielectrics is the subject of in- 
tense on-going interest. While there is no doubt that the 
effect is real, certain suggested applications are somewhat 
controversial. For instance: Schwinger has suggested 
that the Casimir effect might be the underlying physics 
behind sonoluminescence while Carlson, Goldman, 

and Perez-Mercader have suggested possible applications 
to Gamma Ray Bursts More generally, the Casimir 
energy has sometimes been invoked as a possible driving 
mechanism for ultra-high-energy astrophysical processes 
such as quasars. We feel that all aspects of the discus- 
sion could benefit from the improved understanding of 
the basic physics we provide in this paper. 

Historically, the techniques used to investigate the 
Casimir effect were typically a varied mixture of 
Schwinger's source theory, explicit calculations of electro- 
magnetic Green functions (seasoned with time-splitting 
regularization), and sometimes, more physically based 
regulator schemes that take advantage of the analyticity 
properties of the frequency dependent refractive index. 



A key early paper is that by Schwinger, de Raad and 
Milton §. 

Schwinger's most developed point of view can be 
gleaned from the series of papers he recently wrote 
wherein he explored the possible relevance of the Casimir 
effect to sonoluminescence [p]-p[ . For the evolution of his 
views on this subject see [|6|-H]. 

Schwinger found that (for each polarization state) 
the "dielectric energy, relative to the zero energy of the 
vacuum, [is given] by 
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This result can be interpreted in a straightforward man- 
ner as the integral of the difference in dispersion relations 
over the density of states po|JTl[ |. 

In addition to Schwinger's bulk volume term, calcula- 
tions by Milton et al. |[12|-p^ indicate the existence of a 
surface correction. For a dilute (that is, e « 1) spherical 
intrusion of radius R and dielectric constant ei in a di- 
lute dielectric medium of dielectric constant £2 (£2 ~ 1), 
with the eigenmode sum regulated by time-splitting, the 
surface contribution is equivalent to[J 



E 



surface ~ ~T^c(ei — 62)^ -, ■ 

4 [CT)-^ 



(2) 



A controversy has recently arisen over whether or 
not Schwinger's volume term should be retained, and 
whether or not the surface term is the leading term in the 
Casimir energy We have shown elsewhere pO 11 

that the presence of the volume term is generic, and 
have (among other arguments) adduced reasons based on 
density-of-states calculations to bolster Schwinger's cal- 
culation. In this paper we shall pursue this matter fur- 
ther and shall extract as much physics as possible from 
these density-of-states calculations. 

The discussion, though elementary from a technical 
perspective, is quite sufficient to give the most important 



51 ) nf [1^, equation (7.5) of , or the equivalent 
[td. Those calculations only deal with spherical 



^See equation (51 
equation (41) of 

dielectric balls with frequency independent dielectric properties, 
and use an explicit time-splitting regularization. The numerical 
coefficient in this surface term is regularization dependent and it 
does not appear to be possible to relate its absolute normalization 
to the number we will calculate using Schwinger's wave-number 
cutoff. 
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dominant contributions to the Casimir energy. These 
results serve as an important consistency check on more 
sophisticated field-theoretic calculations. 

Furthermore, the present analysis extends Schwingcr's 
result by verifying that generically surface terms do in 
fact show up, but as sub-dominant corrections to the 
dominant volume contribution. General arguments of 
this type are particularly useful because they allow us 
to study arbitrary shapes and not be limited by require- 
ments of spherical symmetry. 

We first discuss some general properties of the bulk 
volume term, noting in particular the dependence upon 
a physically meaningful ultraviolet cutoff, and then turn 
to the issue of finite volume effects. While finite vol- 
ume effects in conductors (or more precisely, for Dirich- 
let, Neumann, and Robin boundary conditions) are well 
understood, the analogous problem for dielectric junc- 
tion conditions (or even acoustic junction conditions) is 
considerably less clear cut. We attack the problem of 
finite volume effects in the presence of junction condi- 
tions via an extension of the Balian-Bloch analysis for 
boundary conditions . We show that the presence 

of a dielectric interface modifies the density of states by a 
term proportional to the surface area of the interface and 
calculate the proportionality constant as an explicit func- 
tion of the dielectric permitivity and permeability. (For 
the related, and simpler, acoustic interface the change in 
density of states is related to the physical fluid densities 
on the two sides of the interface.) 

Finally, we apply this formalism to the estimation 
of the (electromagnetic) Casimir energy in generic di- 
electrics. We show that for dielectric bubbles large com- 
pared to the cutoff wavelength the volume term is dom- 
inant. We point out that the numerical value of the net 
Casimir energy is strongly dependent on the details of 
the high frequency cutoff. Within the context of sonolu- 
minescence this high-frequency sensitivity might explain 
the fact that small admixtures of gas in the bubble un- 
dergoing sonoluminescence can have large effects on the 
total energy radiated: a small resonance in the medium- 
frequency behaviour of the refractive index can be magni- 
fied by phase space effects, and lead to dramatic changes 
in the total energy budget. 

We mention in passing that there will also be an 
acoustic Casimir energy associated with the phonon 
modes The acoustic Casimir energy (while always 

present) is numerically negligible in comparison to the 
electromagnetic effect being suppressed by a factor of 
(speed of sound/speed of hght). 

II. THE DENSITY OF STATES: BULK TERM 

The physics underlying the Casimir effect is that ev- 
ery eigenmode of the photon field has zero point energy 
En = (l/2)?i(jj„; the Casimir energy is the difference in 
zero point energies between any two well defined physical 



situations 

iA\B} = J2-h[uJniA)-^ujn{B)]. (3) 

n 

We always need a regulator to make sense of this en- 
ergy difference, though in many cases of physical interest 
(such as dielectrics) the physics of the problem will au- 
tomatically regulate the difference for us and make the 
results finite. Adding over all eigenmodes is prohibitively 
difficult, so it is in general more productive to replace the 
sum over states by an integral over the density of states. 

Suppose we have a finite volume V of some bulk di- 
electric in which the dispersion relation for photons is 
given by some function coi^k), which describes the pho- 
ton frequency as a function of the wave-number (three- 
momentum) k. Suppose this dielectric to be embedded 
in an infinite background with different dielectric proper- 
ties described by a different dispersion relation uj2{k). We 
regulate infra-red divergences by putting the whole uni- 
verse in a box of finite volume Voo , and calculate the bulk 
contribution to the total zero-point energy of the electro- 
magnetic field by summing the photon energies over all 
momenta (and polarizations), using the usual and ele- 
mentary density of states: [Volume] (i'^fc/(27r)^. (In the 
next section we shall look at finite- volume corrections to 
this density of states.) 

Including photon modes both inside and outside the 
dielectric body the energy of the system is 

/d'^k 1 

/d^k 1 
(^2?^-2(fc). (4) 

Note that outside the dielectric body the photon disper- 
sion relation is that of the embedding dielectric LU2{k). 
Note also that we shall always use the subscript 2 to re- 
fer to the region outside the embedded body, and shall 
use the subscript 1 to refer to the region inside. 

If the embedded body is removed, and the hole simply 
filled in with the embedding medium, we can calculate 
the total zero-point energy as 

f d^k 1 

Ehomogeneous — 2Kx) j 2 '^'^^ ^^-^ 

We define the Casimir energy by subtracting these two 
zero-point energies [00] 

Ecasirnir — E^jjifj^^dQ^—j^Q^y EjiQ^-nogeneous 

/d^k 1 
^-M^i(fc)-^2(fc)]. (6) 

The physical import of this definition is clear: The 
Casimir energy is defined as the change in the zero-point 
energy due to a change in the medium. 
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Note also that the physical meaning of the zero of en- 
ergy is clear: the zero of energy is here taken to be that 
corresponding to a homogeneous dielectric with disper- 
sion relation uj2{k). 

To be obtuse, we could use a different zero for the 
energy — this makes no difference as long as we keep 
the same zero throughout any particular calculation. For 
instance, the zero-point energy of the Minkowski vacuum 
is 



(Pk 1 

E Minkowski = "^^oo / '^27j^ —lick. 



(7) 



Thus an alternative definition for the Casimir energy is 
then 

T^alternative ^ rp 77^ 

Casimir ~ ^embedded— body ^ AHnkow ski 

ffik 1 

W I -h[oji{k)-ck] 



(27r)3 2 



+ 2(Ko - V) 



d^k 1 



(27r)3 2 



h [uj2{k) - ck]. (8) 



For this alternative definition, the zero of energy is clearly 
the Minkowski vacuum. As long as you stick with one 
fixed definition throughout the calculation, or better yet, 
calculate Casimir energy differences directly, quibbling 
about the zero of energy does not matter. (Of course, 
if you change the zero of energy in the middle of the 
calculation the answers will be meaningless.) 

From the general considerations in pO|jll| ] we know 
that the integrand must go to zero at large wave-number, 
and in fact, for any pair of real physical dielectrics the 
integrand must go to zero sufficiently rapidly to make the 
integral converge. 

An integration by parts yields 



Casimir 



vh 

6^ 

Vh 

67r2 



d(fc3) [uj^{k) - UJ2{k)] 



{e [iUi{k) ~ UU2{k)])\ 



[duji{k) - du2{k)]k^ 



(9) 



The boundary term vanishes because of the asymptotic 
behaviour of the uJi{k). The substitution k = u>i(k)ni 
then yields 



Eca 



Vh 



LU-^[n-^{uj) -n-l{uj)]duj. (10) 



While the difference between the refractive indices in the 
above expression goes to zero sufficiently rapidly to make 
the integral converge, it must be noted that the prefac- 
tor of (jj^ implies that the net Casimir energy will be rel- 
atively sensitive to the high frequency behaviour of the 
refractive indices. 



If the Casimir effect ultimately proves to be the correct 
physical explanation for sonoluminescence, this sensitiv- 
ity to the details of the refractive index might plausibly 
explain why sonoluminescence is sensitive to the admix- 
ture of small trace gases into the bubble. (Of course the 
present calculation is static, but the energy calculated in 
this way will be the maximum energy that could possi- 
bly be released in a more realistic dynamical calculation.) 
To make this qualitative statement quantitative we would 
need a detailed model for the refractive index as a func- 
tion of frequency — a task that is beyond the scope of this 
paper. 



III. THE DENSITY OF STATES: 
FINITE- VOLUME EFFECTS 



A. Generalities 

We now look at the contribution arising from sub- 
dominant finite-volume corrections to the density of 
states. The key point here is that the existence of finite- 
volume terms proportional to the surface area of the di- 
electric is apeneric result. The fact that previous calcula- 
tions encountered a surface tension term propor- 
tional to (cutoff)^ is hereby explained on general physical 
grounds without recourse to special function theory. 

We must notice at this stage that the dominant con- 
tribution to the Casimir energy is proportional to vol- 
ume, as the canonical bulk expression for the density of 
states is proportional to the volume: [Volume] d'^k/ {27t)^. 
It is reasonably well-known, though perhaps not so ele- 
mentary, that for fields subject to boundary conditions 
(Dirichlet, Neumann, Robin) the density of states is in 
general modified by finite volume effects. In this paper 
we wish to extend these ideas to fields subject to junction 
conditions (acoustic, dielectric). 

For boundary conditions the general result is 



E 



V 



_fk_ 

(2^ 



s / e 



d^k 
(27r)3/fc 



(11) 



These are the first two terms in an asymptotic expansion 
in l/k. For Dirichlet, Neumann, and Robin boundary 
conditions the coefficients can be related directly to the 
known asymptotic behaviour of the Heat Kernel — they 
are simply the Seeley-DeWitt coefficients in disguise and 
can be obtained, for instance, by suitably transforming 
the results presented in the monograph by Gilkey . 

There are additional terms in this expansion, propor- 
tional to the various monomials appearing in the general 
formulae for the higher Seeley-DeWitt coefficients, but 
we do not further address this issue here except to point 
out that the next term is proportional to the integral of 
the trace of the extrinsic curvature over the boundary. 

An elementary discussion of the general existence of 
such terms can be found in the textbook by Pathria mM , 



3 



while a more extensive treatment can be found in the 
papers by BaUan and Bloch [ibIJi^ . 

For Dirichlet, Neumann, and Robin boundary condi- 
tions the dimensionless variable ^ is a known function of 
the boundary conditions imposed. 

If we let N{k) denote the number of eigenmodes with 
wave-number less than fc, then from the above we can 
write 



N{k) 



1 

2^ 



1 



Vk' + -iSk' + 0[k] 



(12) 



We shall now perform the analogous analysis for junc- 
tion conditions, adapting the Balian-Bloch formalism as 
needed. Our formalism is applicable to both boundary 
conditions and junction conditions. For clarity, and to 
aid in consistency checking, we carry out brief parallel 
computations for the boundary condition case. 



B. Scalar field 

We start for simplicity with a scalar, rather than elec- 
tromagnetic, field. We are interested in the following 
eigenvalue problem 



-I- fc^</. = 0; 



B[4>] = 0. 



(13) 



Here B[(j)\ denotes the boundary conditions imposed. 
Common boundary conditions are tabulated below. 



Dirichlet boundary conditions: 

[(j) — Q on the boundary) 

C = -7r/4. 



(14) 



Neumann boundary conditions: 

{dn4i = on the boundary; where 9„ denotes the normal 
derivative) 



C = +7r/4. 

Robin boundary conditions: 

{dn4i = Kcf) on the boundary; k real) 

C = +7r/4. 



(15) 



(16) 



Surface damped boundary conditions: 

{dncf) = kK(j) on the boundary; k real; note that the eigen- 
value is now explicitly present in the boundary condition 
as well as in the differential equation) 
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1. 



-Im 



In 



1 



1 — in 



arctan(K). (17) 



These results can be read off, for instance, from the 
paper by Balian and Bloch jisj . 

Comparing the Robin and surface damped boundary 
conditions, it might naively be tempting to write 



£,Robin{l^) — S,da7yiped{^/ k) ■ 



(18) 



However in the present context — an asymptotic expan- 
sion in 1/k — such an expression is meaningless. The best 
we can do is to say that 



^Robinin) = lini ^dampedin/k). 



(19) 



Thus for Robin boundary conditions we keep only the 
dominant k ~* oo piece of the Balian-Bloch result. 

On the other hand, in the surface damped boundary 
condition (because of the explicit factor of k appearing 
in this boundary condition) it is meaningful to keep the 
inverse tangent term of the Balian-Bloch result in our 
expression for ^. (As a consistency check, these coeffi- 
cients are also calculated as special cases of the general 
formalism we shall develop below.) 

Acoustic junction conditions: 

We are ultimately interested in junction conditions, 
rather than boundary conditions. For definiteness, we can 
think of an acoustic junction, wherein an acoustic wave 
propagates across some fluid interface: say a bubble of 
some dense fluid embedded in a lighter fluid. In terms 
of the densities of the fluids, (pi, P2), and the velocity 
potentials, ((/)i, (^2), the acoustic junction conditions are 



Pi<Pi 



d„ 



d„ 



(20) 



(21) 



(See |2y, page 24] or gy, page 81]. These two conditions 
represent, respectively, the continuity of the pressure and 
the normal component of the velocity at the interface.) 
We must point out at this stage that the change in prop- 
agation speed and/or density causes a certain amount of 
reflection and refraction, which then changes the density 
of states in the fluid both inside and outside the bub- 
ble {i.e. on both sides of the interface) according to the 
general scheme 



d^k 
(2^ 



s / e 



d^k 

(27r)3fc 



outside 



d^k 

' (27r)3fc 



(22) 
(23) 



For the case of acoustic junction conditions, the dimen- 
sionless variables ^out/in have not yet been calculated. 
We present the calculation below, for now merely quot- 
ing the final result: 



S.out{Pl,P2) = ^ 



Pi - P2 
Pi + P2 



<Upi,P2). (24) 



Formulation of the problem: 

On general grounds, we expect the ^ to be a function of 
both the acoustic refractive index (that is, a function of 
the relative acoustic velocities), and the relative densities. 
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If vq is some arbitrary reference speed, we can define the 
refractive indices by 



ni = kivo/uj 



and 



n2 = k2VQ/uj, 



(25) 



and further define the relative refractive index by n = 
ni/n2. (Note in particular that w is continuous across 
the interface, whereas ki is not.) It is also useful to define 
the density contrast by p = pi/p2- 

In the special case where there is no dispersion, the 
phase and group velocities are equal and we simply have 



ni 



vq/vi and n2 = Vo/v2- 
We know, from first principles, that as n 



(26) 



1 and 

p — > 1 the acoustic boundary becomes indistinguishable, 
as both fluids have the same density and refractive index, 
so we must have 



£.out/tn{n,p) 



as 



1 and p 1. (27) 



To calculate ^(n, p) for acoustic junction conditions, we 
modify the discussion of Balian and Bloch page 407] 
to derive an expression for ^(n, p) in terms of an integral 
involving the reflection coefficient R{p, n; k). 

We start from the result for the density of states in 
terms of the time-independent Green function [|l5|, equa- 
tion (II. 6), page 409]. Taking N{k) to be the number of 
modes with wave-number less than fc, we can construct 
a suitably smoothed density of states formally described 
by the relation 



Pdos{k) = 



dN 



dk 



(28) 



smoothed 



(Details of the smoothing procedure can be found in p5| .) 
Note that we prefer to express the density of states in 
terms of the wave-number ki rather than in terms of the 
variable E — kf. See equation (1.3) on page 402 of [ p5[ . 
Thus 



dN dEdN , „r,, , 

In terms of the asymptotic expansion of interest 
Pdos{k)^^{Vk^ + S^k + 0[l]). 



(29) 



(30) 



Working on either side of the interface (with i taking on 
the values "in" or "out" as appropriate) equation (II. 6) 
on page 409 of (l5| yields 



Pdosi^i) 



2k- f 

— d^x \imJm[G{x,x';h + ie)], (31) 

TT J x' — *x 



where the integration over x now runs only over region i 
as appropriate. 

It is important to realize that the Balian-Bloch formal- 
ism is built up under the assumption that all the eigen- 
values are real — this constrains the type of problems we 
can deal with to loss-free undamped situations. 



We are interested in an arbitrary interface, but pro- 
vided the interface is smooth, we can locally replace it 
by its tangent plane. This approximation is equivalent to 
neglecting sub-dominant pieces proportional to the trace 
of the extrinsic curvature. (If we were interested in ex- 
plicitly calculating the next coefficient in the expansion 
we would have to locally approximate the surface by its 
osculating ellipsoid, as done for the case of boundary con- 
ditions by Balian and Bloch.) 

Truncating the expansion at the surface area term, we 
locally approximate the interface by a plane interface, 
located at z = 0, with region 2 (the outside) at z > 
and region 1 (the inside) at z < 0. We are only interested 
in the diagonal part of the Green function. To calculate 
this diagonal part in region i we can assume the source is 
also in region i and write the total Green function in this 
region as a sum of a direct and a reflected contribution. 

The direct part of the Green function is responsible 
for the bulk contribution to the density of states, while 
the reflected part of the Green function gives the surface 
contribution. Since, in the tangent plane approximation, 
we are dealing with a perfectly flat interface higher order 
contributions are explicitly excluded. 

The volume contribution has already been calculated 
in and we are now interested in the extra piece 

of the Green function that arises from reflection at the 
interface. Using cylindrical coordinates, the contribution 
to the Green function due to the reflected wave can be 
put into the Sommerfeld representation (an integral over 
transverse wave- number kt) 



— / R\h,kt)Mktr) 
47r Jo 

exp[iK{ki,kt){z + z')] 



(32) 



(See equation (4.2.5) on page 103 of |^^, with an ap- 
propriate change of notation.) Note that R(ki, kt) is the 
reflection coefficient. It is a function of the frequency 
and the transverse wave-number and will consequently 
depend on the precise nature of the boundary conditions 
imposed. The Sommerfeld representation has the inter- 
esting feature that it expresses a Green function, which 
is related to the behaviour of spherical waves, in terms 
of a reflection coefficient defined for plane waves. Here 



Kih,kt) = 
More explicitly 

Kout{k2, h) 
Kin{ki, kt) = 



kl 



kl 



kl 



kl 



(33) 



(34) 



(35) 
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If we look at the diagonal part of this reflection contribu- 
tion, (x = x'), and note that Jo(0) = 1 we immediately 
see 



K'ih, kt) eyi^[2iK{ki, kt)z\ 



in 



Kih,kt) 



ktdkt. (36) 



(Note that we are calculating what is in field theory par- 
lance an off-shell Green function. The integration over 
kt is an integration over all off-shell transverse momenta, 
and this integration is not to be limited by any on-shell 
constraint such as kt < fc^.) 

For the density of states (counting only the appropriate 
contribution arising from either side of the interface, that 
is, z > or z < 0) 



r re flection 

ik. 



S" / dz Im 
27r^ Jo 



kt dkt R^{ki + ze, kt) 



exp[2i_R'(fci A- ie, kt)z] 
K{ki+ie,kt) 



(37) 



The z integration is trivial. (Because ki has a small pos- 
itive imaginary part, which is inherited by K[ki,kt), we 
can guarantee convergence of this integral.) 



Preflectioni^i) 



47r2 



S Im 



W{h+ie,kt) 
K{ki + ie, ktY 



kt dkt 



(38) 



It is useful to define the dimensionless variable u = 
kt/(ki -\- ie), so that u has a small negative imaginary 
part. We get 



Preflectioni^i) — ^^"^ 



W{ki, u — ie) 
1 - (w - »e)2 



du 



(39) 



If we now take this contribution to the quantity pdosi and 
convert to the ^' variable as defined in this paper using 



we find 



Equivalently 



-Im 



RHk,+ie,kt) 
K{ki + ie,kt)^ 

°° R'{ki,u-ie) 
1 - (u - «e)2 



kt dkt 



udu. 



(40) 



(41) 



(42) 



This is our general result for the surface contribution to 
the density of states. The surface term is seen to be a 



suitable average of the reflection coefficient appropriate 
to the boundary conditions at hand. (Note that if we 
were on-shell, we would interpret u as the sine of the 
angle of incidence, and u would then be limited to the 
range u S [0, 1]. As this is an off-shell computation for 
the off-shell Green function, the range of integration goes 
all the way to infinity and trying to interpret u as the 
sine of the angle of incidence only leads to unnecessary 
confusion. Indeed, in calculating this Green function, we 
are effectively dealing with a spherical incident wave, so 
there are many angles of incidence 9i. To identify u as 
the sine of the angle of incidence only makes sense for 
an incident plane wave, and is in the present context 
meaningless.) 

The application of this result to specific cases of inter- 
est merely requires us to calculate the relevant reflection 
coefficients and perform the integrations. 

The integral for standard boundary conditions: 

In some well known cases the relevant integrations are 
straightforward. For example for Dirichlet, Neumann, 
and Robin boundary conditions the reflection coefficients 
are —1, -f 1, and -1-1 respectively, and integrating out to 
some large cutoff value of u we have 



1 



1 - (u - ie)"^ 



udu 



1 



\ — {x — ie) 



dx 



^-[-\n{l-{x~^e)}]\", 

= i (ln{l + ie) - ln{l - + ie}) 
= -[i7r + ln(C/2 _ iy^i2 



(43) 



Note that the integral itself diverges, though the imag- 
inary part is both finite and independent of the cutoff. 
Taking this imaginary part gives 



(44) 



This reproduces the standard results quoted above. 
[Equations (0-|l6|).] 

The surface damped boundary condition is a little 
trickier. In this case the reflection coefficient can be 
shown to be 



R{u) = 



(45) 



(See, for example, equations (3.4.4) and (3.4.5) on page 
87 of DeSanto |2^, and translate to our notation. Note 
that an analytic continuation in k is required to turn the 
surface impedance boundary condition discussed there 
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into the surface damped boundary condition discussed 
here.) 

Subtracting and adding 1 to the integrand converts the 
integral into that encountered in the previous calculation 
plus an integral that is well-behaved at infinity. The 
relevant integral is again elementary: 



1 



\/l — + in 



udu = +2 In (1 + in) . 

(46) 



Taking the imaginary part of the above reproduces the 
result announced in equation (pT[): 



^ = — — arctan(K). 



(47) 



Checking the above: 



l-{u-itY 



ier — IK 



udu 



-1 

oo 



1 




(u2- 


-le) 


— IK 


1 - 




~ ie 


) 




(m2 - 


-^e) 


■f iK 




1 










-2iK 




1 - 


- (w2 - 


- ie 


) 




(u2- 




-f iK 




1 








iK 




du' 


1 - 


- u' + 


ie 


Wl-u' 


+ 


ie 


+ iK 




1 






in 






du" 




u" 


— ie 


w- 


-u" + ie 


+- in 






1 






K 




du" 




u" 


— ie 


Wu" 


~ ie + K 





udu 



udu 



u — le 
1 



K 



u + K — te 
1 



du 



In 



u — le u + K — le 
u — ie 



du 



= -2 



In 



u + K — te 
—i — ie 



-I + K — te 
= +21n(l + iK). 



(48) 



(The original contour was chosen to run underneath 
the two branch cuts emanating from u = —l + ie and 
u = +1 + ie; thus under the change of variables u" = 
vz?— 1 the branch cut must be chosen so that the new 
contour terminates at —i and not at +i.) 

The integral for acoustic junction conditions: 

We are finally ready to study the case of interest: acoustic 
junction conditions. The reflection coefficient is now pO| , 
equation (3.f.l9), page 82] 



R{p, n; u) 



py/l — u^ 



(49) 



Consistency check I: Note that p —t +oo gives i? = +1, 
as appropriate for Neumann and Robin boundary condi- 
tions; p ^ gives i? = — 1 as appropriate for the Dirich- 
let boundary condition; while p —f oo with k — —in/p 
fixed gives the surface damped boundary condition. 

Consistency check II: Similarly n — > -f-oo gives R = 
— 1, as appropriate for Dirichlet boundary conditions; fi- 
nally n — > +ioo gives R = +1 as appropriate for Neu- 
mann and Robin boundary conditions. 

Observation: The reflection coefficient exhibits an in- 
version symmetry as we move from one side of the inter- 
face to the other, this symmetry being inherited by the 



Rinip, n; u) = i?o«t(l//3, 1/n; u). 



(50) 



Thus 



Cm(p, n) = Co„f(l/p, 1/n). (51) 
Calculation: We are interested in evaluating 
udu p\J\ — u^ — \J"n?- — u^ 



— Im 



1 



p\J\ — u^ + V "n? — u^ 



(52) 



The integrand has a pole at w = 1 of residue —1/2, and 
branch cuts emanating from u = ±1 which can be chosen 
to terminate at u = ±n. Asymptotically, as m ^ oo, the 
integrand goes as 



lp-1 



u p 



1 



(53) 



This is already enough to tell us that the imaginary part 
of this integral can be finite if and only if p is real. For the 
acoustic equations this is actually very sensible physically 
since it is meaningless to drive the density complex. To 
evaluate this expression we subtract and add 1 to the 
integrand, and make use of the integral J udu/{l — u^), 
evaluated in equation (E3[), to write 



TT 

+ 2 
+ Im 



udu I pvl 



1 -u2 



(54) 



This conveniently gets rid of the pole so that the integral 
is now unambiguously finite. Indeed 



Q = 



TT 

2 



2plm 



udu 



vr 



(55) 



Now we also have to take n to be real, otherwise we step 
outside the Balian-Bloch formalism. For now, also take 
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n > 1, the alternative case being completely analogous. 
The integrand is now imaginary only over the range u € 
[1, n], and we can change variables to set 



Q 



That is 



du' 
l-u' 



-Im 



\/n^ — u' + ip\Ju' — 1 



Q 



du' 



1 I ^ u' {v? - u') + p^{u' - I) 



Equivalently 



Q 



n -1 J„," 



du" - 1 - u"V^ 



Now define u" = {n^ — l)w. 



Q 



^ dw y/1 — wJw 



(56) 



(57) 



(58) 



(59) 



Note that the refractive index n has now completely dis- 
appeared from the integral. This gives 



0- 



2 



P- 



1 



1 



1 



(60) 



We can re-do the calculation for n < 1. A few interme- 
diate steps change but the final result is the same. We 
finally have our announced result 



(61) 



TT 


' P - 


1' 


TT 


Pi 


- P2 


4 


.P + 


1 


~ 4 


.Pi 


+ P2_ 



Note the remarkable result that this is independent of n 
for n real. With hindsight, we can see that the acous- 
tic junction conditions explicitly make reference only to 
the density of the fluid, and not to the velocity of sound 
(refractive index), which might be viewed as an a pos- 
teriori justification for the absence of refractive index in 
the final result. However we know of no simple physics 
argument that would justify this, and must rely on the 
explicit calculation presented above. 

As p ^ +00 we recover Neumann and Robin boundary 
conditions while as p —^ we recover Dirichlet boundary 
conditions. Also note that on interchanging the two re- 
gions, p ^ so we have 

^tn{p) = £,out{l/p) = ~(,out{p), (62) 
as expected from our earlier discussion [see Eq. (pl|)]. 



C. Electromagnetic field 

For the electromagnetic field, we can use the analysis 
presented by Balian and Bloch in pages 273-274] to 
view the electromagnetic eigenvalue problem as a combi- 
nation of vector and scalar eigenvalue problems. A stan- 
dard result is 

Perfect conductor boundary conditions: 

{E X n ^ and i? • n = on the boundary) 



(63) 



This vanishing of the surface term for perfect conduc- 
tor boundary conditions is due to a cancellation between 
TE and TM modes. (For a surface of general shape the 
separation into TE and TM modes is meaningless; TE 
and TM modes make sense only in situations of extremely 
high symmetry. Nevertheless, sufficiently close to any 
conducting surface we may approximate the surface by its 
tangent plane — and in this approximation the decompo- 
sition into TE and TM modes makes sense. The general 
vector minus scalar decomposition alluded to above then 
approximately reduces to the simpler scalar plus scalar 
decomposition for the TE and TM modes.) 

Dielectric junction conditions: 

For the case of ultimate interest we are of course inter- 
ested in dielectric junction conditions. A full apprecia- 
tion of the (perhaps unexpected) subtleties involved with 
dielectric junction conditions might be gleaned from the 
fact that even for a plane interface the situation is suffi- 
ciently complicated as to warrant a recent 600 page tech- 
nical monograph [ pl| , and a continuing stream of research 
papers (see for instance @). 

Nevertheless we can make a few general statements on 
physical grounds before doing a detailed calculation of 
^. In analogy with the case of the scalar field, finite- 
volume effects will distort the density of states both in- 
side and outside the dielectric body according to the gen- 
eral scheme 



E 

inside 

E 

outside 



V 



d^k 

(2^ 



s / e 



{Voo - V) 



(2^ 



d^k 
(27r)3fc 



d^k 
(27r)3fc 



(64) 
(65) 



For the case of a dielectric junction, we expect ^(e, p.) to 
be a function of the permeability and permitivity, and 
we know, from first principles, that as e —^ 1 and p —^ 1 
the dielectric boundary disappears as both media become 
the same, so we must have 



C(e,p)->0 



1 and 



1. 



(66) 



When we turn to including dispersive effects we note that 
^(e, p) should ultimately be taken to be a function of the 
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wave-number dependent quantities e(fc), fi{k). Since we 
know that as — > oo the dielectric must ultimately mimic 
individual atoms embedded in vacuum, we must have 



as 



(67) 



The calculation of ^ for the electromagnetic field is 
an easy exercise given our results for the acoustic prob- 
lem. We decompose the electromagnetic field near the 
approximately plane boundary into TE and TM modes. 
In terms of the relative refractive index, relative permi- 
tivity, and relative permeability, the reflection coefficients 
(for the outside region) are simply ||2^, pages 83-84] (or 
see equations (86.4) and (86.6), page 295], or p^ , 
pages 281-282])! 



=(p = fi,n;u) 



(68) 



B fields is sufficient to completely determine the electro- 
magnetic field. In terms of these normal components the 
junction conditions are simply 
TE mode: 



TM mode: 



dnH" — dnH2- 



eiEi — e2E. 



2 ) 



(73) 



(74) 



(75) 



(76) 



Applying the formalism derived for the acoustic junc- 
tion conditions, the previously quoted results for ^ im- 
mediately follow. 



eVT 



+ ^^r 



(69) 



(Remember that n = ^JT\i. Also, we have defined n = 
ni/n2, e = ei/e2, and /i = ^1/^2) 

Thus, applying the previous acoustic results, we get 
the remarkably simple formulae: 



/-TE/ 







TT 


>i 


- M2" 


4 


/i + 1. 


^ 4 


.Ml 


+ M2. 


TT 


"e- 1" 


TT 


ei 


- £2 


4 


.£ + 1_ 


^ 4 


.ei 


+ £2. 



C'^ip)- (70) 



^IrTie). (71) 



Note that the result for the TE modes is independent of 
e, while that for the TM mode is independent of /z. Since 
most typical dielectric materials are magnetically inert, 
/i « 1, the TE contribution is typically much smaller 
than the TM contribution. 

Consistency check: Instead of appealing to the iden- 
tification of reflection coefficients, we can get the same 
results directly from the dielectric boundary conditions. 
We know that 



E' 



eE"-' 



H^, and fiH' 



(72) 



must be continuous across the boundary. 

If we are dealing with a plane interface, or in the ap- 
proximation that we are sufficiently close to a curved in- 
terface, specifying the normal components of the E and 



Be careful with all the different notations in use. 



IV. THE CASIMIR ENERGY 

Including these surface contributions to the density of 
states, the total zero-point energy for a dielectric body 
embedded in a background dielectric is easily seen to be 



E. 



embedded— body 



= 2V 



T^o^^l(^) 



(27r)3 2 



25 / ^hc 



(27r)3 2 
+ 2{Voo - V) 

+ 2S 



6ri(e,A*) 



^i(fc) 
ck 



:-huJ2[k) 



(27r)3 2 



d^k 1^ 
: — nc 



(27r)3 2 



£,out{e-,y) 



UJ2{k) 

ck 



(77) 



This is just the generalization of equation (|j) above to 
include surface effects. The quantity ^ denotes an average 
over TE and TM modes. To calculate the Casimir energy 
we now simply subtract the homogeneous dielectric zero- 
point energy [equation (Q)] to obtain 



Eca 



= 2V 



2S 



d^k 1 



(27r) 
d^k 1 
(27r)3 2 



3 ^H^iik) - uj2{k)] 



he 



Cm(e,/^) , ?o«t(e,/^) 



Til 



"2 



(78) 



Even though the surface terms seem to be additive, there 
is a "hidden" minus sign, as we shall see below, due to 
the fact that Cm(p) = -S.out{p)- 

This is quite enough to give a good qualitative feel for 
the physics: the Casimir effect will in general induce a 
surface tension that goes as (cutoff )3. 
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It is useful to define 



S(ei,/Lfi;e2,Ai2) 



ni 



n2 



(79) 



and so write the Casimir surface tension as 



CT (surface tension) 



(2^ 



he S(ei,/^i; £2,^2)- (80) 



From our previous results for ^, taking the case of mag- 
netically inert media for simplicity (/i = 1), we see 



^(ni,n2) = - 



1 

ni 



1 

f^2 



(81) 



Here we indeed see that the two surface terms contribute 
with opposite signs, largely canceUing each other. We 
can factorize this to yield 



^(711, 712) = 



TT {rii - n2Y{ni + 7^2) 
8 nin2(n\+n\) 



(82) 



Note that this vanishes as (rti — 772)^, with one factor 
of {til — 712) coming from the fact that the individually 
tend to zero as ni — > rt2 and the second coming from the 
partial cancellation discussed above. 
What does this do to the Casimir energy? 



Casijnir 



= 2V 



2S 



(fk 1 , 
(27r)3 2 

d^k 1 ^ TT (rji 



n2 - ni 



77.177,2 



"2)^('^1 + "2) 



(27r)3 2 8 nin2{nl+nl) 



(83) 



This is our general result for the Casimir energy. We 
now insert a momentum dependent refractive index into 
the above to explictly evaluate the coefficients. The phys- 
ical cutoff is provided by the fact that both refractive 
indices are known to tend to 1 at large momenta. 

A naive hard cutoff, following the ideas of Schwinger, 
simplifies these expressions considerably. Naive hard cut- 
offs are of course an idealization that suppresses much of 
the physical detail, and are justified only for order of 
magnitude estimates and for comparison with the previ- 
ous literature where naive hard cutoffs are often the only 
extant results. Suppose we take 



and 



ni(fc) = m &{K -k) + e{k - K), 



n2{k) ^ 772 Q{K -k) + Q{k - K). 



(84) 



(85) 



(It is an additional gross over-simplification to set the 
cutoffs for the two media equal to one another, but it is 
standard and is the only way to make connection with 
previous calculations. Keeping separate cutoffs for the 



two media is straightforward but algebraically somewhat 
messy.) 

The Casimir energy is then given by 



1 

8^2 
1 



Ecasimir = ^ 



1 1 

TLl 772 

£.in{ni,n2) , S,out{ni,n2) 



ni 



77-2 



while the Casimir surface tension is 



-he K-^ 



6n("-l,"2) , £.out{ni,n2) 



771 



«2 



(86) 



(87) 



Inserting the specific formulae for ^ then yields 



AStt 771772(77,2 -I- 772) 



Now particularize to dilute media, by taking 771 
1 



E, 



Casimir 



8tt 
1 

487r 



^VheK^ [772 



771 



S he [(771 - ^2)2] -I- 



772. 



(89) 



The volume term here is the dilute medium limit of 
Schwinger's result Q, while the surface area term repro- 
duces the Milton et al. result JT^-p^ . There is an overall 
normalization difference between this surface term and 
the special case calculated by Milton et. ai, this normal- 
ization difference being attributable to a different choice 
of regulator. The critical physics lies in the volume ver- 
sus surface area dependence, the power of the cutoff de- 
pendence, and the behaviour as a function of refractive 
index. 

Note that the bulk term is dominant if 

V K:$>S, (90) 
that is, for dielectrics with linear dimensions satisfying 



L^{V/S) > 1/X = Ao/(27r). 



(91) 



For a typical dielectric we estimate Aq « 1 , 000 
Angstrom, so for dielectrics of this size or greater the 
Casimir energy will be dominated by bulk effect. This 
is eertainly the ease for sonoluminescenee where typical 
bubble radii are of order 100, 000 Angstroms. 

For small enough dielectric particles the surface term 
will not be negligible in comparison to the volume term — 
this is no great surprise to people studying mesoscopic 
systems for which the existence of finite volume effects is 
well known. 

Finally, we mention that for [non-dispersive] Neumann, 
Dirichlet, and Robin boundary conditions the existence 
of a surface term contributing to the total Casimir energy 
has been known for some time — see for instance [Eql. 
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V. DISCUSSION 

The main results of this paper are: 
(1) The Casimir energy in a dielectric medium is domi- 
nated by a volume term. Indeed, for a finite-volume of 
dielectric 1, embedded in an infinite volume of different 
dielectric 2, 

f (Pk 1 

EcasiniiT = / 77:-^-?i [wi(fc) - i^2{k)\ 



2S 



(27r)3 2 
d^k 1 
(27r)3 2 



?ic^(ei,^i;e2,Ai2) 



(92) 



where the dots represent terms arising from higher-order 
distortions of the density of states due to finite-volume 
effects. 

(2) If we adopt a simple cutoff model for the dispersion 
relation, the volume term is 



1 

ni 



1 

n2 



(93) 



This result is completely in agreement with Schwinger's 
calculation in |]l| , and in disagreement with ||l^,|ljl . 
(3) In addition, there will be a sub-dominant contribution 
to the Casimir energy that is proportional to the surface 
area of the dielectric. This surface contribution takes the 
generic form 



E, 



surface 

Casirnir 



1 



She 



(ni - n2)^(ni + n2) 



nin2{n{ + n^) 



(94) 



This term is sub-dominant provided 
l//5> l/i<r = Ao/(27r). 



(95) 



(4) In general, we can expect these to be the first two 
terms of a more general expansion that includes terms 
proportional to various geometrical invariants of the 
body. By analogy with the situation for non-dispersive 
Dirichlet, Neumann, and Robin boundary conditions p^ ] 
we expect the next term to be proportional to the trace 
of the extrinsic curvature integrated over the surface of 
the body. 

(5) The analysis of the present paper has been limited 
to situations of real refractive index (loss-free insulating 
dielectrics). Generalizing to lossy conducting media is 
clearly of interest but will require a careful re-assessment 
of the entire formalism. 
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